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Abstract 



A simplified theoretical model for the constant-magnitude steady-state surface 
photovoltage (SPY) technique is developed. Emphasis is placed on the determination of the 
minority carrier diffusion length L in the particular case for which the sample thickness is 
quite lower than the true value of L. The model is derived from the theoretical basis with the 
assumption that incidental light flux varies linearly in the explored range. This model fits very 
well to experimental data points obtained from both FZ and CZ mono-silicon crystals. We 
showed that in the case of thin wafers, SPY technique is not adapted to L measurements. It 
leads only to extracting a short length proportional to the sample thickness. However, if the 
recombination velocities and the minority carrier diffusion coefficient are known, then we 
could reach the true value of L by a computer calculation. 
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I. INTRODUCTION 

The surface photovoltage (SPY) technique is a well-established contact-less 
experimental tool for the characterisation of semiconductors, which relies on analysing 
illumination-induced changes in the surface voltage. Recently, a single treatise devoted solely 
to an in-depth description of this technique has been published [1]. It provides a tutorial 
review covering nearly all aspects of the technique, from the earliest theoretical and 
experimental achievements to the latest developments. 

One of the most important and successful applications of the SPY technique is the 
determination of the minority carrier diffusion length L in the bulk. This application has been 
introduced by GOODMAN [2,3], and well developed by both MOORE [4] and CHIANG et 
al. [5]. It is certainly the branch of SPY measurements, which has had the most significant 
impact on the semiconductor manufacturing industry. 

Using SPY to determine L parameter has numerous advantages and large success over 
others techniques [1]. However, it is not frequently employed because of many quantitative 
limitations of the approach and the validity of assumptions used [1]. One of the necessary 
experimental conditions required by SPV-based diffusion length measurement states that the 
sample thickness w must be much larger than L parameter. PHILLIPS [6] has shown that this 
condition should be met for samples thicker than four times their L. Otherwise, the extracted 
length is only an effective value, which is not trivially related to any real diffusion length. This 
effective value should be shorter than the true one because of the recombination of minority 
charge carriers at the backside of the sample [6]. For much thinner samples, the extracted 
effective length tends toward a limit value, which should be quite lower than the true value of 
L. From experiments undertaken by EICHHAMMER [7], this limit value was found to be not 



2 



much larger than half the sample thickness. But, no explanations related to this particular 
result have been given. This work is mainly devoted to focus on this point and to try 
conducting a deep examination of such approaches of the problem. In this aim, we propose a 
simplified theoretical model and experiment on the relationship between the sample thickness 
and the extracted minority carrier diffusion length by SPY technique. 

II. MODELLING 

Let us consider an n-t)^e semiconductor sample with a w thick and a large front- 
surface barrier. When monochromatic light, with wavelength A, and an incident photon flux 
density 4>o, is directed into the surface barrier, the semiconductor absorbs it provided that the 
photon energy is near or above the band-gap energy of the material. As a result of photon 
absorption, electron-hole pairs are created. The optical generation rate G(z) of these electron- 
hole pairs at a distance z from the illuminated surface is [8]: 

G(z) = aOo(l-R)xexp(-az) (1) 

where a is the optical absorption coefficient of the semiconductor, R is the optical reflection 
coefficient at the illuminated surface and % is the quantum photon efficiency for electron-hole 
pair generation. Some of the generated electron-hole pairs recombine immediately at surface 
edges of the sample. To take this into account, a lower photon quantum efficiency T| 
substitutes for % in equation (1). Also, while diffusing, a part of the generated charge carriers 
recombine through complex mechanisms which we can integrate in a simple exponential 
function as exp (-z /L) where L is the minority carrier diffusion length in the bulk of the 
sample. 

Based on the above, the hole flux and, then, the photo-courant density Jp resulting 
from the photo-generated pairs in the whole w thick of the sample may be expressed as: 
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q a Oo(l-R) Ti e-«2 e-2^L (2) 



where Ci is a constant and q is the absolute value of the electron charge. 

In general, a, R and r| all are functions of A,. Within the wavelength range of interest, 
however, R and r| are practically constant. The proportional term (l-R)Ti in (2) will always 
appear together with 4>o. For simplicity, we therefore write ^>o(1-R)t| as ^. So, (2) becomes: 

Jp = Ci Tq a a> e-(«+VL)z (3) 

which, after integration, gives: 

Jp = CiqO^^[l-e-(i+«L)w/L] 
1 + aL 

In the well known first case [2] assuming that the w » L condition is satisfied, the 
equation (4) is reduced to: 

Jp = CiqO-^ (5) 
1 + aL 

As the open-circuit photovoItageVoc developed at the illuminated surface barrier is 
given by [9]: 

nkT^ fj ^ 



Voc = Log 

q 

then, we can write: 



P+1 



(6) 



Jp = Js t'l^-/"''^ - l] = CiqO^ 



1 + aL 

where Js is the saturation current (Js ~ qDpo/L), n is the ideal factor, which is generally 
close to 1, po is the bulk hole density at thermal equilibrium and D is the hole diffusion 
coefficient. 
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Keeping Voc at a constant value by adjusting and hence 4> while varying A, and thus 



a(k), we obtain: 



O = C2 



Dpo 



L 



aL 



= C2 



Dpo 



L + 



a 



(8) 



where C2 is also a constant: C2 = Ci (ei^"^/"'''^ - l) 



(8') 



This equation is similar to that (expression number 30) formulated by CHIANG [5] 
when the term wa is very small. From this equation (8), L can be determined by extrapolating 
the 4> versus 1/a plot to ^> = 0; the negative of the intercept with the 1/a axis yields L. Of 
course, it is assumed that a(k) is known over the required range. 



From the same equation (8), we can also easily draw a convenient expression for L as 



follows: 



L = 



Oiai - <l>2a2 O1X2 - O2X1 



(<t)2 - Oi) ai ai (O2 - Oi) 



where: xi,2 = l/cci,2 (9) 



Let us consider now the second case for which w « L. From our knowledge, this 
particular case has never been well-discussed or/and developed. To consider it, the principle 
of the method consists in assuming that the flux variation, which maintains the phototension 
Voc at a constant value, is linear in the exploration zone allowing writing: 



rn 



= o 



f 1 ^ 



d$ 



+ 



a 



a ai 



(10) 



or: 



0(x) = <D(xi) + 



dO(x) 



dx 



(x-xi) where x= 1/a (10') 



Xl 
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By extrapolating the 4>(x) plot to ^> = 0, the negative of the intercept with the x axis yields i 

0(xi) 



I = 



dO(x) 




dx 


Xl 



Xl 



(11) 



with 4>(x) can be expressed from the equation (4) as follows: 

, ^ ^ Jp X + L 1 
0(x) = — 



L 



w w 

1 - e "^7 



(12) 



It should be noted that the first term of the second member in the equation (11) is, in 
fact, the inverse of the logarithmic derivative of 4>(x) at xi: 



0(xi 


) 




dO(x) 




d 


dx 


Xl 


dx 



LogO(x) 



Xl 



XI + L 



+ 



w w 
- + — 
w e I L Xl 



w w 

1 - e I L Xl 



(13) 



When we put this in the equation (5), it becomes after deduction 

_ x^L - (xi w + wL + Xl l) Xj e 



£ = 



w w 

— + — 
L Xl 



w w 

x2 + (xi w + wL - x^ )e [l Xl 



(14) 



Let us perform now a development in series of the exponential function up to the 
second order as follows: 



w ^ w 

e 1 L X, I ^ 1 - 



rw 
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Xl 


J 



(15) 



When we put this in the equation (14), we accomplish a simple but a little long and tiresome 
process to lead to: 



(16) 
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This result, remarkable by its simplicity, means that in any case where the sample 
thickness is quite smaller than the minority charge carrier diffusion length L, the measured 
lengths i using SPY technique are nothing of the kind. They only give an adequate value of 
the sample thickness. 

It could be interesting to mention that CHIANG formulas [5] are less different to those 
we are proposing in this work. This difference should be due to the more explicit details given 
by this author on the minority carrier losses by recombination processes. Indeed, from our 
side, these losses have been entirely integrated in a simplified expression which leads not only 
to the same conclusions on the determination of L, but also specify the sense to give to the 
SPV measurements when using thin wafers. 

III. EXPERIMENTAL VALIDATION 

To valid the remarkable theoretical result expressed by equation (16), some 
experimental measurements have been performed in the required corresponding conditions. 
They have been conducted on virgin, n-tj^e, phosphorus doped, <100>, 1-5 Ilcm mono- 
silicon crystals from Wacker. Both Czochralski pulled (CZ) and float-zone grown (FZ) single 
crystals have been used. They present a very similar initial true diffusion length L with a 
predetermined average value around 500 |jm. Samples of dimensions typically 1x1 cm area 
were cut out from circular wafers of 3 in. diameter and 450 fim thick. These samples were 
first lapped on both faces with Aloxite 125 grit (average grain size = 9 |jm) to perform serial 
deeper and deeper sectioning allowing to obtain a series of thickness w ranging from 150 to 
nearly 400 |am in variable steps as measured with a micrometer. Next, one face (i.e., the front 
side) of each sample was slightly etched in (5 : 3 : 3) 70% HNO3 : 50% HF : 100% 
CH3COOH chemical mixture for 60 s, which would remove roughly 15 |im thick. After, all 
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prepared samples were degreased in trichlorethylene, cleaned in acetone, carefully rinsed in 
running de-ionised water and finally dried under a nitrogen gas flow to be ready for the SPY 
measurements in the particular w « L experimental condition leading to extract I effective 
values instead of L. 

i was measured on the etched face of the sample by means of a home built SPY set-up 
using a small pickup probe (0 = 1 mm). To have a good SPY signal, all samples were boiled 
in de-ionised water for 1 hour to ensure correct band bending at their front surface. The 
relative error in extracted I value is estimated at 10 %. Just before the measurements, the 
samples were cleaned again to minimise external contamination and remove the native oxide 
layer. 

Figures 1 and 2 show the dependence of extracted I on the sample thickness in the 
case of FZ and CZ crystals respectively. Each I value reported in these figures has been 
obtained from repeated measure performed on nearly the same area of the sample. The straight 
line fit to FZ data points and that of CZ, show very similar slopes. From these experimental 
results, it is very easy to point out that, at least in the chosen thickness range, the measured 
length i varies quite proportionately with the sample thickness w according to that expected 
and predicted by the theoretical model described above. As shown, we find slopes around an 
average value of 0,5 for the two kinds of crystals. This result agrees very well with that 
explain-less empirically deduced from the experiments undertaken by EICHHAMMER [7]. 

Using the whole equation (14), when we plot the evolution of I as function of the 
sample thickness w (figure 3), one can show clearly that for some L, w and xi values, 
corresponding roughly to our experimental conditions, the linear approximation used in the 
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model is well justified. We find then slopes of 0,5 average value, with small variations to be 
associated, in less explicit way, with the recombination mechanisms and the true diffusion 
lengths. As we can see in the figure 3, within the framework of our assumption, we should 
have to use samples with thickness at least three times larger than the diffusion length to hope 
performing a correct measurement of this parameter by SPY technique. 

Although the curve slope depends partly on the diffusion length value, we saw that its 
variations remain small and could not lead to any determination of L except to determine, in 
other respects, parameters such as D, po as well as recombination velocities in the bulk and at 
the involved surfaces of the sample. 

IV. CONCLUSION 

In this work, another dimension in the influence of the wafer thickness on the minority 
carrier diffusion length L measurement by SPV technique emerges. It is shown that in the case 
of wafer with low thickness compared to L, SPV technique was not adapted to L 
measurements since it amounts to extract a length proportional to the sample thickness. If, by 
another way, the values of the recombination velocities in the space charge region as well as 
on the involved surfaces of the sample and D coefficient are known, then L should be 
obtained by using computer processing either by CHIANG model [5] or from equation 14. 
This requires first measuring the optical reflection and absorption coefficients, the doping 
density in the sample, the open-circuit photovoltage at the illuminated surface and the incident 
photon flux density. Moreover, it needs to make simple assumptions on the low-level 
excitation and on the linearity of 4> function. 
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FIGURE CAPTIONS 



Figure 1: 

Dependence of the SPV extracted / value on the FZ mono-silicon crystal wafer thickness. The 
straight line fit to data points shows a slope of about 0,57. 

Figure 2: 

Dependence of the SPV extracted £ value on the CZ mono-silicon crystal wafer thickness. The 
straight line fit to data points shows a slope of about 0,55. 

Figure 3: 

Dependence of the calculated /value on the semiconductor wafer thickness. For some L, w 
and xi values, corresponding roughly to our experimental conditions, the linear 
approximation used in the model is well justified. We find then slopes of 0,5 average value. 
We should have to use samples with thickness at least three times larger than the diffusion 
length to hope performing a correct L measurement by SPV. 
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